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Abstract 

Polymetric walls are walls built from bricks in more than one size. 
Architects and builders want to built polymetric walls that satisfy cer- 
tain structural and aesthetical constraints. In a recent paper by de Jong, 
Vinduska, Hans and Post these problems are solved by integer program- 
ming techniques, which can be very time consuming for patterns consisting 
of more than 40 bricks. Here we give an extremely fast method, generating 
patterns of arbitrary size. 
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1 Introduction 

How to built a brick wall with different sizes of bricks? This problem was 
recently discussed in the paper |2]. 




Fi gure 1: An example of a polymetric brick wall from www.metselen.net/fotoalbum.html 

Let us try to formulate the requirements put forward by architects and 
builders. These are 



1 



[Rl] The absence of regularity in the patterns, as e.g., periodicity 
[R2] Long horizontal joints, to facilitate the building process 
[R3] The absence of long vertical joints, for sturdiness and strength 
[R4] Pleasing aesthetics, as e.g., not too many large bricks close together 

One can appreciate the fulfilment of all four requirements in Figure [l] 

In [5] the four requirements are met by considering the construction of a 
brick wall pattern as a Pallet Loading Problem, which is then solved with a 
Mixed Integer Linear Program. Here we take a completely different approach, 
constructing the patterns as pseudo-self-similar tilings generated by 2-D substi- 
tutions. Well known properties of such tilings automatically yield the fulfilment 
of [Rl], and also more or less [R4]. Requirement [R2] can very easily be built 
in, so that [R3] becomes the main challenge. 

The style of this note will be somewhat informal, for precise definitions and 
more context see the papers [5J and [5J. See Section 14 (and the book cover) 
of Allouche and Shallit's monograph |T] for the subclass of tilings generated by 
multi-dimensional automatic sequences. Since then the literature on tilings has 
expanded a lot (cf. [6]). 

2 Prouhet-Thue-Morse patterns 

The Prouhet-Thue-Morse sequence 0,1,1,0,1,0,0,1,... is the one-dimensional se- 
quence generated by the substitution — !• 1, 1—^10. It is folklore that in the 
"same way" one generates a 2D-sequence (see, e.g., [3]) by the 2D-substitution: 

10 1 

i ^ 1 , 1 h-s- 1 0. 

Usually this sequence is then represented in the plane by white and black 
squares. When we take 2x1 and 3x1 bricks this corresponds to the following 
substitution rules 




The third iteration of this substitution starting from the 2x1 brick B21 
yields the following pattern 
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Obviously requirement [R3] will be violated, since we obtain vertical seg- 
ments crossing the whole pattern. To cope with this, we slightly change the 
previous substitution to the substitution r given by 




Now a fourth order iteration of r starting from the 2x1 brick yields the pattern 
in Figure [2] 




Figure 2: Skewed 2D-Prouhet-Thue-Morse: r 4 (B 2 i). 



These skewed P-T-M patterns nicely satisfy [R3], see the following proposi- 
tion. 

Proposition 1 For any n the vertical segments in the patterns T n (B2i) and 
T™(I?3i) have length at most 2. 

Proof: By induction on n, using the following figure. 




□ 



3 General patterns 

We would like the pattern to be generated by a geometric 2D-substitution, 
mapping rectangles i x j (bricks By) to unions of rectangles (bricks), resulting 
in a self-similar tiling (see [6] for the basic definitions). However, as already 
clear from the first P-T-M example above, such self-similar tilings will violate 
the bounded vertical segments condition [R3]. Thus we should rather con- 
sider pseudo-self-similar or pseudo-self-affine tilings as studied by Priebe and 
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Solomyak ([!]). In a self-affine (rectangular) tiling with expansion coefficients 
Ai and A2, all the bricks (usually called proto-tiles in the literature) are scaled 
by a factor Ai in the a;-direction and a factor A2 in the y-direction, and then are 
exactly tiled by translations of the original set of bricks. Note by the way that 
the P-T-M-substitution is self-similar with two lxl bricks (expansions Ai = 2 
and A2 = 2) , but not with a set consisting of the 2x1 and a 3 x 1 brick. 

In a pseudo-self-affine tiling the bricks do not tile exactly the expanded 
bricks, but have bricks 'sticking out' and 'sticking in', in a way that still preserves 
the tiling property. We will consider a subclass of 2D-substitutions which have 
the desired properties, which we call of PP-type, i.e., which have the periodic 
parallelogram property. With this we mean that the image of a i X 1 brick is a 
region as in Figure [3J left. The black horizontal segments have length iAi, and 
the two curves on the left and the right are translated copies of each other. For 
a i x j brick this is repeated periodically j times, cf. Figure [3j right. 



Figure 3: Left: Image region for a i x 1 brick. Right: Image region for a i x 2 
brick. 

A simple example of a PP-type substitution with Ai = A2 = 2 and the three 
brick alphabet A = {£?n, -B21, -B22} is given by 




The patterns a 3 (Bn), a 3 (B 2 i) and a 3 (B 22 ) are 




Let v max (ij,n) be the maximal length of a vertical segment in a n (Bij) for 
n > 1 and B t j e A. For the substitution above v max (ll 7 1) = 1, v max (21, 1) = 
2,f m ax(22, 1) = 3, and one can prove that for this substitution 

«max := sup{w max («j,n) : B i3 € A, n > 1} = w max (22,3) = 11, 
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Figure 4: a 4 (B 2 2). 



but in general it is not an easy task to determine the largest possible length 
u m „ of a vertical segment (which, of course, can be infinite). 

The following approach generates patterns satisfying [R3] in an easy way. 
We say the pattern a(Bij) has no crossings if there are no vertical segments 
consisting of boundaries of the bricks in a(Bij) that connect two points on the 
boundary of the region determined by a(Bij). In the substitution above, for 
example, a(Bn) and o(Bi2) do not have a crossing, but a(Bi\) does. 

Proposition 2 Let a be a W-type substitution with expansions Ai and A2. 
Suppose that for all Bij the pattern a(Bij) has no crossings. Let j* be the 
maximal height of the bricks. Then any vertical segment in any pattern a n (Bij) 
has length at most 2j*(\2 — 1). 

Proof: Let [(x, y), (x,y')] be a vertical segment. The points (x,y) and (x,y') 
lie in some a(Bij) and a(Bi'j'). If these are the same, then \y — y'\ < j*(\2 — 1). 
If these are different, then the no crossing property implies that they are neigh- 
bors of each other, and the length of the segment is at most twice the maximal 
height minus one of a a(Bij), i.e., at most 2j*(A 2 — 1). □ 

We give an example with Ai = 2 and A2 = 3. 



■ ^ 

This substitution is of PP-type, and has no crossings. According to Propo- 
sition^ (with j* — 1), the length of vertical segments in any pattern generated 
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by this substitution is at most 4. Actually a refinement of the proof of Propo- 
sition [2] shows that is even bounded by 3, and this bound is achieved as can be 
verified in Figure [5] 




Figure 5: A clip from ct 5 (_B 2 i). 



We end this section with a few remarks on the number of bricks that are 
needed in a pattern. To this end one introduces the matrix M a , with entries 
M a (ij, i'f) which by definition are equal to the number of times the brick B^y 
occurs in the pattern cr(_By), for all B^ and B^y e A (in the literature some- 
times the transpose of M a is taken as definition). 

Proposition 3 The maximal eigenvalue of M a is A1A2. The asymptotic fre- 
quencies of the bricks are given by the normalized left eigenvector corresponding 
to this eigenvalue. 

Proof: This proposition is well-known, see e.g. [5]. Here we give a short 
proof, in a slightly more general situation. We will assume that just as in the 
self-similar case a satisfies for all B^ G A 

Area(cr(Sy-)) = A1A2 Area(-Bij). 

Note that our PP-type substitutions satisfy this equation. Now let v be the vec- 
tor with entries Area(-Bij), and w the vector with entries AresL(a(Bij)). Then, 
since Area(cr(i?ij)) is the sum of the areas of the Byy times the number of times 
Bi'ji occurs in a(Bij), we obtain 

M a v = w = X1X2V. 
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Conclusion: v is a right eigenvector of M a with the eigenvalue A1A2. But since 
all the entries of v are strictly positive, A1A2 has to be the Perron- Frobenius 
eigenvalue of M a . The second statement in the proposition is a well-known 
result (see e.g. [7]), generalized by Peyriere (see [5]) to higher dimensions: the 
normalized left eigenvector of M a corresponding to the Perron-Frobenius eigen- 
value equals the asymptotic frequency of the bricks in large patterns. □ 

As a final remark we mention that the PP-property can be extended in the 
obvious way to substitutions that in addition have a periodic parallelogram 
structure in the horizontal direction. We illustrate this with an example: 




The expansions are here both 4, there are no crossings, and it is easily shown 
that f m ax = 3. The patterns <7 2 (-B2i) and cr 2 (i?22) are given by 




l__J l_ l ' U U 1 




Li Ll 1 l_ l l__J ' 



4 Random brick wall patterns 

The interpretation of the requirements [Rl] (absence of regularity) and [R4] 
(pleasing aesthetics) is of course for a large part a matter of taste, and some 
might consider the pattern in Figure [5] too 'rigid'. One way of making the 
patterns less rigid is to introduce randomness in the construction. As an in- 
troduction we start with a very simple self-similar example, which hence does 
not satisfy requirement [R3]. We take the 1x2 and 2x2 bricks with random 
substitution rules given by 



I 
I 



Here random means that each expanded brick is replaced by one of the two 
possibilities with probability 1/2, independently of the other replacements. In 
this simple example the substitution matrix M a is non-random, since the two 
possibilities are just permutations of each other. A simple computation learns 
that for n = 1, 2, . . . 

2 . 4"- 1 4™- 1 



1 4^ 2 • 4 n— i 

Note that different sequences of coin flips will lead to different patterns, so there 
will be 

21+6H _ 22-4 n -l 

different realizations of a n+1 (-B22) , since the total number of bricks in cr"(i?22) 
is equal to the sum of the two entries in the lower row of M™. This means that 
the three realizations in Figure^are picked from a set of 2 2 ' 43-1 = 2 127 « 10 39 
elements! 

■ ■ ■ ■■ ■ ■ ■■ ■ ■ ■ ■ ■■ ■■ ■ ■ ■ 

■ ■■■■■■■■■ ■ ■ ■ AS" 



■ ■■ ■ 



■■■ ■ ■ ■ ■ ■ ■■■ ■ ■ 



Fi gure 6l Three realizations of cr^ (B22) • 



In general the matrix M a will be random, and its expectation will give us 
the mean growth of the number of bricks in the patterns a n (Bij). The theory 
of multi-type branching processes can be applied to obtain further information 
on these patterns. 

We next consider a pseudo-self-similar random substitution, where the coin 
is flipped with a success probability p for some p € [0,1]: 

■ 

The maximal length of a vertical segment now is a random variable V max , 
and in this pseudo-self-similar example the behavior of V max (p) as a function of 
p is particularly interesting. Note that Vm a x(0) = 00, but that V max (l) = 2. 
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Figure [7] shows 2 realizations for p = 1/3 and p = 2/3 each. 
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